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I.  INTRODUCTION 

Vertical  sound  velocity  gradients  in  the  ocean  are  several 
orders  of  magnitude  larger  than  horizontal  sound  velocity 
gradients.  Upon  this  "proposition",  horizontal  refraction  of 
sound  propagation  is  neglected  in  most  underwater  sound  ray 
tracing  models.  However,  in  regions  of  sharp  horizontal 
thermohaline  gradients,  such  as  oceanic  fronts  and  eddies, 
horizontal  refraction  may  not  always  be  disregarded  for  long- 
range  sound  propagation.  This  note  investigates  the  importance 
of  horizontal  refraction  through  cold  and  warm  eddies. 

Water  masses  such  as  the  Sargasso  Sea  or  Slope  Water  provide 
a simple  setting  in  which  to  study  sound  propagation  because  of 
their  large  areal  extent  and  weak  horizontal  gradients.  Only 
vertical  density  stratification  of  these  water  masses  affects 
refraction  patterns  of  sound  ray  paths  (neglecting  bottom  effects) . 

Recent  papers  by  Gotthardt  (1973)  and  Parker  (1971)  however, 
reveal  that  Gulf  Stream  rings  (eddies)  of  different  thermohaline 
properties  may  be  common  throughout  the  northwestern  protion  of 
the  Sargasso  Sea  and  in  the  Slope  region  between  the  Grand  Banks 
and  Cape  Hatteras.  The  physical  properties  of  these  rings 
(typically  200  km  in  diameter)  differ  from  the  physical  properties 
of  the  water  mass  surrounding  them.  Cold  eddies  have  been  shown 
to  alter  the  vertical  refraction  sound  ray  paths,  (Gemmill  and 
Khedouri,  1974,  Vastano  and  Owens,  1973),  but  little  attempt  has 
been  made  to  understand  whether  an  oceanic  eddy  will  alter  sound 
propagation  in  horizontal  plane  (a  problem  frequently  referred 
to  as  bearing  error) . Horizontal  sound  speed  changes  across  these 


NAVOCEANO  TN  6150-24-74 


eddies  have  been  calculated  to  be  20-30  m sec-1  in  50  km 

(Khedouri  and  Gemmill,  1974).  Although  these  sound  velocity  j 

changes  represent  small  gradients,  the  relative  refractive  index* 

of  the  eddy  and  surrounding  water  is  about  0.98  for  a warm  eddy 

and  1.02  for  a cold  eddy  (between  air  and  water  the  relative 

index  of  refraction  is  1.33,  i.e.).  These  values  are  nearly 

unity  but,  at  long  ranges  horizontal  refraction  of  sound 

propagation  may  be  important. 

A model  to  examine  the  horizontal  refraction  of  sound  rays 
through  a cylindrical  eddy  using  geometric  principles  of  optics 

was  developed.  In  a physical  sense,  the  model  can  be  interpreted 

■ . } 

to  portray  the  amount  of  horizontal  refraction  caused  by  an  eddy 

on  sound  propagation  between  the  ocean  surface  and  the  deep  .1 

\ '1 

sound  channel.  The  model  results  are  compared  with  the  results 
of  the  thin-lens  approximation.  Keep  ii  mind  that  this  model 
isolates  only  the  horizontal  plane,  and  in  no  way  minimizes  the 
importance  of  vertical  refraction.  | 

II.  THE  MODEL  AND  HORIZONTAL  REFRACTION  PATHS 

This  model  was  developed  assuming  an  idealized  cylindrical 
eddy.  The  effect  of  eddy  currents  on  sound  ray  paths  are 
assumed  to  have  at  least  an  order  of  magnitude  less  than  the 
effect  of  horizontal  sound  speed  differences.  The  geometry  of 

the  cylindrical  eddy  is  shown  in  figure  1.  | 


*relative  refraction 

water  is  defined  as  VQ/VE  where  V, 
surrounding  water  and  Ve  is  the 


index  of  an  eddy  and  its  surrounding 
0 is  the  speed  of  sound  in  the 
speed  of  sound  in  the  eddy. 
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The  center  of  an  eddy  of  radius  (R)  is  at  a distance  (Xo) 
from  a point  sound  source  located  at  the  origin  (o)  x=o,  y=*o. 

Rays  are  drawn  at  1°  intervals  between  -40  and  +40°  at  the 
source  (where  0 is  measured  positive  counterclockwise  from  the 
x-axis  to  the  ray) . Each  ray  that  intersects  the  eddy  is 
refracted  upon  entering  the  eddy  and  on  leaving  the  eddy  to 
some  predescribed  distance  behind  the  eddy.  Given  the  sound 
speed  outside  (VQ)  and  inside  (Vg)  the  eddy,  the  refracted  ray 
paths  are  determined  from  Snell's  law.  The  derivation  is  presented 
in  Appendix  A. 

Horizontal  refraction  diagrams  were  plotted  for  both  cold 
and  warm  eddies.  Figures  2 and  3 show  the  refraction  patterns 
through  a cold  eddy,  figures  4 and  5 through  a warm  eddy  with 
radii  of  100  km;  for  sound  sources  placed  at  200  km  and  at  an 
infinite  distances  from  the  eddy  center.  The  sound  speed  is 
set  at  1500  m sec”1  in  both  eddies.  Sargasso  sound  speed  was 
given  a value  of  1530  m sec  1 and  Slope  Water  sound  speed  1470 
m sec"1  (relative  refractive  index  of  1.02  and  0.98,  respectively). 
By  setting  the  sound  source  at  infinity  the  sound  rays  will 
converge  to  (diverge  from)  a region  called  the  focal  length  of 
the  cold  (warm)  eddy. 

Cold  eddies  produce  two  important  alterations  to  normal  ray 
paths.  Refraction  through  the  cold  eddy  reduces  spreading  of 
the  rays,  but  more  important,  shadow  zones  are  formed  down- 
range  from  the  eddy  edge  (figures  2 and  3) . 

-4- 
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A warm  eddy  acts  to  disperse  sound  rays  more  rapidly  as 
they  refract  through  the  eddy  (figures  4 and  5) . It  is  possible 
that  not  all  rays  that  intersect  the  warm  eddy  will  refract. 

Rays  intersecting  the  eddy  near  its  edge  may  have  an  angle  of 
incidence  greater  than  the  critical  angle,  and  will  be  reflected. 

III.  COMPARISON  TO  THIN  LENS  OPTICAL  MODEL 

The  spatial  variation  of  the  image  produced  by  an  eddy  was 
examined  from  combinations  of  the  eddy  radii  (R)  and  the  eddy 
acoustic  contrast*  (AC)  with  sources  at  varying  distances  from 
the  eddy  center  (Xo) . These  computations  were  compared  to  the 
results  produced  using  the  thin— lens  equation  (Jenkins  and  White, 

* W .>  r -7  •'?.  •*.♦•-{  I ' ,9 

1957),  and  some  of  the  comparisons  are  presented  in  Table  1. 

1 X 2AC 

The  thin-lens  equation  can  be  written:  (— -+  — ! > 

where  S is  the  distance  of  the  source  to  the  lens,  and  S1  is 
the  distance  of  the  image  from  the  lens.  Note,  for  negative 
distances,  the  images  "appear"  to  be  in  front  of  the  eddy.  For 
most  combinations  of  eddy  radii,  and  acoustic  contrast  and  source 
distances,  images  will  appear  to  be  between  the  source  and  the 
eddy.  For  images  to  occur  behind  eddies,  the  source  would  have 
to  be  outside  the  focal  length,  or  on  the  order  of  1000  km  away 
from  the  eddy.  The  comparison  between  the  thin-lens  image 
calculation  and  the  cylindrical  eddy  image  model  calculation 
indicates  that  the  thin-lens  equation  is  a good  approximation. 


♦acoustic  contrast  is  defined  as  (Vo”Ve)/Ve  where  VQ  ir 
the  sound  velocity  outside  the  eddy  and  Ve  is  the  sound  velocity 
within  the  eddy. 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

Sound  rays  will  spread  horizontally  more  rapidly  upon 
refraction  across  a warm  eddy  compared  to  totally  homogeneous 
water  mass.  The  reverse  is  true  for  a cold  eddy,  resulting  in 

a shadow  zones  down  range  from  the  eddy  edge. 

Focal  lengths  and  position  of  apparent  images  are  approximately 

the  same  using  either  this  cylindrical  model  or  the  physical  thin- 
lens  equation. 

It  appears  that  large  horizontal  gradients  in  the  thermohaline 
structure  of  the  ocean  should  not  be  neglected  during  sound 
propagate  .dies.  Field  experiments  to  determine  the  importance 

of  hori  -fraction  are  planned. 


KH'I  (Cold  Eddy  with 


Figure  2 Eddy  Radi 


TABLE  I 

COMPARISON  OF  CYLINDRICAL  EDDY  IMAGE  AND  THE  THIN-LENS  IMAGE  POSITION 
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APPENDIX  A 

THE  HORIZONTAL  REFRACTION  MODEL 


The  eddy  horizontal  refraction  model  was  developed  assuming 
an  idealized  cylindrical  eddy.  Sound  speed  is  specified  as  a 
constant  within  and  outside  the  eddy.  The  geometry  for  the 
model  was  shown  in  Figure  1. 

A point  sound  source  is  located  at  the  origin  (o)  x=o.  y-o, 
which  is  at  a distance  (Xo)  from  the  center  of  an  eddy  of  radius 
(R) . Sound  speed  outside  the  eddy  is  given  as  VQf  within  it  is 
VE.  For  a given  initial  angle  (8  measured  positive  counter- 
clockwise from  the  x - axis) , the  point  of  intersection  for  a 
ray  from  a source  at  origin  to  the  eddy  can  be  found  by  solving 
simulataneous  equations  of  a straight  line  (from  the  source  to 
the  circle)  and  of  a circle  (the  eddy  of  radius  R)  at  distance 
Xo  from  the  source.  The  two  equations  are: 

Y = X tan  8 equation  of  straight  line  (la) 

(X-Xo)2  + Y2  = R2  equation  of  circle  (lb) 

By  substituting  la  into  lb  to  eleminate  Y,  and  rearranging 
the  terms  and  returning  to  la,  the  sound  ray  will  intersect  the 
eddy  at: 

2 


Xj  = xo  cos2e+  xo  cos4e  - (x3  - r2)  cos2e 
Yj  = Xj  tan  6 


(2a) 

(2b) 
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At  the  point  of  intersection  the  ray  will  be  refracted  by 
the  eddy  according  to  Snell's  law: 


sma 


where  VQ  and  oQ  are  the  sound  velocity  and  the  incident  angle 

at  intersection  with  the  eddy,  and  the  E subscripted  terms  are 

the  values  after  refraction  in  the  eddy.  The  incident  angle 

(aQ)  is  the  angle  between  the  incident  ray  and  the  normal  to 

the  circle.  The  angle  of  intersection  to  the  eddy  is  simply 

computed  from  the  slope  of  the  normal,  n 

Yt 

of  the  intersecting  ray,  mR  = 
relation: 


Xj-Xq  • ana  SJ-U1 

(4b) , using  the  trigonometric 


The  refracted  angle  through  the  eddy  is  computed  rearranging 
Snell's  law 


The  ray  will  be  refracted  again  upon  leaving  the  eddy.  The 
coordinates  at  the  point  of  ray  termination  relative  to  the 
point  of  exit  from  the  eddy  are  computed  from: 

AXt  = D cos  (0  - 2A0)  (10a) 

AYt  * D sin  (6  - 2A0 ) (10b) 

or  the  termination  point  relative  to  the  origin. 

X,,  = Xj  + AXe  + AXj  (11a) 

Yt  = Yj  + AYe  + AYt  (lib) 

The  set  of  equations  above  can  be  reduced  for  the  special 
case  of  a source  placed  an  infinite  distance  away.  In  this  case, 
the  source  acts  as  an  infinite  "line  source"  so  that  all  the  rays 
are  assumed  parallel  to  the  x - axis  and  can  be  placed  any 
distance  from  the  eddy.  The  governing  equation  can  be  written  as: 

(X-Xc)2  + Yb2=R2  (12) 

For  this  case  XQ  is  the  arbitrary  distance  from  the  "line  source" 
to  the  vertical  axis  through  the  eddy's  center  and  Ys  is  the 
perpendicular  distance  to  the  x - axis.  Since  Ys  = Yj  (parallel 
ray  to  x - axis)  then  Xj  is  solved  from  (12)  using: 

Xj  - p0  1 r2  " ^/2  (13a) 

Yj  - Ys  (13b) 

.-15- 
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Obviously  (13)  does  not  hold  if  Yg  >_  because  rays  do  not 
intersect  the  eddy  but  pass  straight  on  by.  For  Ys  <R, 
equations  (13)  are  used  in  place  of  equations  (2)  to  find  the 
point  of  intersection  of  the  incoming  ray  with  the  eddy.  The 
computations  following  the  ray  paths  through  the  eddy  are 
calculated  from  (3-12)  as  before. 
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